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D . Abstract: In this paper, we consider the problem of quantum measurement of neutrino 



velocity. We show, that the well known neutrino flavor oscillations are always accompanied 
by the oscillations of neutrino velocity. In particular, the velocity of a freely moving 
neutrino is demonstrated to periodically exceed the speed of light. Unfortunately, the 
superluminal effect turns out to be too small to be experimentally detected. It is also 
shown that neutrino velocity significantly depends on the energy, size and shape of the 
neutrino wave packet. Owing to the big experimental error of the recent experiments, 
these dependences remained unnoticeable. Finally, we have shown that the recent claims 
that superluminal neutrinos should loose energy during their flight is not true. Instead, 
our formula suggests the approximate conservation of energy along neutrino trajectory. All 
these results have been obtained without violation of special theory of relativity. 
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1. Introduction 

The existence of neutrino flavor oscillations is considered a well established fact in con- 
temporary physics |2|, ^, ^, ^, ^ ^, |9[. It is caused by the fact that the three known 
neutrino flavors, fg, i'^ and Vr, are not the eigenstates of the Hamiltonian. Instead, they 
are superpositions of the true eigenstates vi, V2 and z^s, having sharp masses mi, m2 and 
777-3, respectively. As a consequence, the flavors v^., and Ur oscillate during the time 
evolution of a free neutrino. 

The purpose of this work is to examine if similar oscillatory character may appear 
when it comes to neutrino velocity. It is motivated by the observation that flavor oscilla- 
tions necessarily imply the oscillations of neutrino masses. Then, owing to the momentum 
conservation, we expect a freely moving neutrino to have oscillating speed. 

The idea that neutrino velocity is closely related to the flavor oscillations is not new. 



It has been explored in refs lid, 11, M], with the result that neutrinos are superluminal. 



However, our analysis shows that the main expression of refs |1C, 11, f2l represents just 



a small correction to our result. In particular, we demonstrate that neutrino velocity 
considerably depends on the size and shape of the neutrino wave packet. 

In this paper, we shall work in the approximation of just two flavors: and Vr- We 
shall demonstrate that the velocity of the free muon neutrino indeed has oscillating char- 
acter. In particular, the neutrino velocity periodically exceeds the speed of light. The 
maximal value of the detected velocity along neutrino trajectory has also been found. In- 
terestingly enough, the probability to detect maximal velocity is pretty small, and in some 
cases, goes to zero. At the same time, the obtained formulae turn out to be unexpectedly 
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sensitive to the shape and size of the neutrino wave packet. This makes the comparison 
with the known experiments very difficult. For one thing, the exact shape of the experimen- 
tal wave packets is not known. For the other, we are not convinced that all the packets in 
the ensemble are identical. Nevertheless, we have tested our formula by comparing its pre- 
dictions with three recent experiments |13, For that purpose, the numerical values 



of our free parameters (such as neutrino energy) are chosen from these experiments, and 
the undetermined free parameters, related to the wave packet shape and size, were chosen 
by consulting the literature [^ ]. As a result, a good agreement with related measurements 
has been achieved. In particular, the derived energy dependence of the neutrino velocity 

been shown to remain undetectable in the considered exper- 
iments. The same holds for the apparent loss of energy during the flight of super luminal 



neutrinos |17|. Indeed, we have shown that the rate at which superluminal neutrinos loose 
their energy is linear in time, but the slope of the graph E(t) is extremely small. This way, 
the loss of neutrino energy becomes unnoticeable in all terrestrial experiments. 
In what follows, we shall use the natural units h = c = 1. 

2. Neutrino dynamics 

To simplify the study of neutrino oscillations, in what follows, we shall adopt two useful 
approximations. The first is that only two flavors, and i^r, will be considered to span 
the internal Hilbert space of the neutrino (the remaining internal degrees of freedom will 
be neglected). The second is that the neutrino moves along the rr-axes, and its dependence 
on y and z is considered irrelevant. In practice, this means that we reduce our task to a 
one-dimensional problem. With this in mind, the Hilbert space of the neutrino becomes 
T-L = C"^ Ho, where is the internal Hilbert space spanned by two orthonormal vectors 
(i = 1,2), and "Hq is orbital Hilbert space spanned by the momentum eigenvectors 
\p), (— oo < p < oo). The basis vectors \i)\p) are taken to be the eigenvectors of the 
Hamiltonian of the free neutrino: 



H\t)\p) = Ei\i)\p), Ei = ^mj+pK (2.1) 

Thus, the states \i)\p) have the sharp values of mass. On the contrary, the muon and 
tau neutrinos are not the eigenstates of the Hamiltonian. Precisely, the general and v-^ 
states have the form 

= j dpa{p){cose\l)\p) -sine\2)\p)) , (2.2) 
\vr) = j dpb{p){sme\l)\p) +cose\2)\p)) , (2.3) 

where 6 is the mixing angle determined by 0.92 < sin^ 29 < 1, and 

J dp\aip)\^ = J dp\b{p)\ 
ensures the proper normalization. 
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In what follows, we shall consider the temporal evolution of initially pure muon neu- 
trino. Its generic state vector is given by (|2.2[ ), and its evolution is determined by the 
Hamiltonian (|2.1|). Thus, we obtain 



\v^{t)) = j dpa{p)(^cosee-'^^^\l)\p) -sinee-'^^^\2)\p) 

The probability density to detect the muon neutrino in point x is then given by 

P^{x,t) = \{u^{x)\u^{t))\\ 

where \v^{x)) = cosO \^)\x) — sm.6 \2)\x) stands for the eigenstate of the position operator 
for muon neutrino. Direct calculation then yields 



j dp a{pyP'^ ( cos^ e e-*^i* + sin^ 6 e'^^^t 



(2.4) 



To solve this integral, we shall use the reasonable assumption that the distribution a{p) is 
sharply localized around the value p = Po- This way, the neutrino is characterized by an 



almost sharp value of momentum. Then, we can expand the energy Ei{p) = ymf 

in a power series around p = po, and keep only linear terms. We shall also introduce the 

group velocities 

dEio 

Vi = , 

dpo 

and the shorthand notation EiQ = Ei{pQ). With these approximations, the integration in 
dU) leads to 

Pf,{x, t) = Ai cos^ 9 e-i(^^o-viPo)t + ^2 sin2 6 e-*(^2o-^2Po)< ^ ^ (2.5) 



where Ai = A{vit — x), and 



A(t) = , / dpa(p)e 

V2^ J 



IpT 



The amplitude A{t) in the coordinate space is the exact Fourier transform of the momen- 
tum amplitude a{p). We have already assumed that the momentum of the neutrino wave 
packet is well localized around p = pq. To take this explicitly into account, we shall take 
the neutrino wave packet in the form 

Air) = /,(r)e— P» , (2.6) 

where the modulus p{t) is localized around r = 0, and the wavelength 2it/pq is much smaller 
than the packet size. The latter ensures a small uncertainty of the packet momentum. In 
what follows, we shall see how our results depend on the size and shape of the neutrino 
wave packet. 

The values of energy and velocity of the two terms in (|2.5| ) differ as a consequence of 
different masses they carry. To estimate their difference, we introduce 

_ £^20 - Eio A 2 _ 2 2 

oj — , At; ^ V2 — vi , Am = m2 — rrii . 
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In the ultrarelativistic limit po ^ rrii, suitable for the description of neutrinos, Eio is further 
decomposed as 

Eio = po + -^ + --- io = —— + ■■■ , 2.7 
and the velocities Vi take the form 

mf ^ 2uj , , 

Vi = l--A2+--- => Av = + • • • . 2.8 

2po Po 



Let us now analyze the probability density ( |2.5D in more detail. First, observe that the 
velocities vi and V2, although close to the speed of light, are different from each other. 
Owing to this, the two initially overlapping packets will gradually separate. After a long 
enough time, we shall see two distinct neutrino wave packets. To estimate the time needed 
for the separation of the two packets, we make use of the packet size. Then, the time 
needed for their minimal separation (when the distance between the packets reaches the 
packet size) is given by 

2e _ Mpl 



t 



\Av\ Am? 



where £ is half the size of the wave packet. Using the numerical data from the recent 
experiments |13, 14, ^] (as shown in Table || of the last section), we find 



t > 5 



mm, 

telling us that the two wave packets in (|2.5D practically coincide in any terrestrial experi- 
ment (5 min ~ 10* km). For this reason, in what follows, we shall simplify our considera- 
tions by adopting the restriction 

t\Av\ <^i. 

Note that this restriction still allows for very long flights {x ~ 10^ km). 

Let us evaluate the amplitudes in the first order in the small parameter tAv. To this 
end, we introduce a new time coordinate r, and the average velocity v, defined by 

+ ^ , - _Vl+V2 . . 

t=- + T, v= — - — . (2.9) 
V 2 

The new time coordinate measures time relative to the moment a particle with the average 
velocity v arrives at x. Having in mind that v ~ 1, we can say that r measures the neutrino 
delay as compared to the arrival time of the photon. The moduli of the amplitudes now 
become: 

Pin t-x)= p{vT - —t) pivT) yP (^'^) ' 

/ \ /- Av . ,_ . tAv ,,_ , 

p{V2 t-x)= p{vT + —t) W p{vt) + —p [vt) . 

Owing to the smallness of the factor | | , the second term is expected to be much smaller 
than the first one. Indeed, the numerical value of this factor in the recent experiments 
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P^ , |T5| is less than 10 At the same time, the phase of A, as defined in ( |2.6D , is 
subject to no approximation at all. 

With the adopted approximations, the probability density ( |2.5| ) takes the form 

Pf,{x, t) = p^{vt) (1 - sin^ 29 sin^ ujt) - [p^{vt)]' cos 26 , (2.10) 

where vt = vt — x, and the prime denotes derivative with respect to the argument. As the 
realistic wave packets have finite size, it is natural to assume that the amplitude /?(r) is 
localized in the interval —i<T<i. This way, the time coordinate r is restricted by the 
packet size, which is typically much smaller than x. In what follows, we shall adopt the 
reasonable restriction 

\T\<i<.X. (2.11) 



3. Velocity oscillations 

In this section, we shall study the motion of muon neutrino by studying the spacetime de- 
pendence of its probability distribution (p.lOf ). To this end, let us first place the neutrino 
detector in a fixed position x. This way, the probability to detect the muon neutrino be- 
comes a function of time, only. The moment neutrino arrives at the detector is determined 
as the time the probability density ( 2.1(\) reaches its maximum in the point x. The needed 



time of arrival is then obtained by solving the equation 
With the approximation v ~ 1, it ultimately leads to 



p2 p2 • p 

with 

Lof sin^ 26 sin 2u}t 



T2{t) + -Tl{t)=0, (3.1) 



4 l-sin2 20sin2a;t + ^cos20' ^^'^^ 

T2(t) = -- t. (3.3) 

Pol- sin^ 26 sin^ ojt + ^ cos 26 



Before we continue, note that the restriction ( 2.11 ) allows us to easily switch between the 



t and X coordinates. Indeed, after taking into account |r| <C x, and v ^ 1, the equation 
(U) shows that 

t ^ X . 



With this, every solution of ( |3.lD , which is originally in the form r = T{t), can approxi- 
mately be rewritten in the form r = t(x). 
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3.1 Example of a simple wave packet 

The above analysis is the most we can do without specifying the amplitude A. As it turns 
out, the solution r = t(x) depends a great deal on the size and shape of the wave packet. 
Before we carefully examine this dependence in the next section, let us describe the simple 
example of a particular wave packet. Up to the normalization constant, we define the 
modulus of the amplitude (|2.6|) as 



p(t) oc < 



, otherwise , 



where the parameter £ determines the packet size. To simplify calculations, we shall work 
in the approximation t <^ I. First, we multiply the equation ( |3.1| ) with £^ — r^, and 
expand it in a power series of r/£. Then, we drop terms proportional to t^/£^, whereupon 
the equation (|3.lD becomes the quadratic equation 



(n + 2t2) + £V - (n + T2)t^ = . (3.5) 
As seen from (3.5), Ti — )• implies r — )• 0, which uniquely determines the physical solution 



of the quadratic equation to be 



On the other hand, if we restrict our considerations to the linear approximation in r/£, the 



equation (|3.5|) yields 

T{x)=Ti{x)+T2{x). (3.7) 

What is immediately seen is that ti(x) is a periodic function of x with the period 

L = -. (3.8) 

(For illustration, the neutrino energy of 17 GeV yields L ~ 18,000km). The term T2{x) 
is linear in x, and the coefficient in front of it is periodic with the same period L. In fact, 
this term is the main result of Refs. [0, |T^. It does not depend either on the size or 
on the shape of the wave packet. Owing to the linear dependence on x, the term T2 can 
exceed ri in the limit x — )• oo. Let us accurately estimate the conditions needed to prevent 
T2 from dominating ti. Thus, we start with \t2\ *C |ti|, and obtain 

„ Pn^^ I sin 2ujx\ 

T - ^ « ^ ■ 3.9 

4 8 X 

Generally, this condition can not be fulfilled for all distances. However, if we avoid distances 
which are too close to nL/2, the requirement ( ^.9]) can be solved for x. For demonstration 
purposes, let us take data from the recent experiments |13, 14, 15], as shown in Table |2|. 



Then, one finds ^ 1.15 • 10^ km, and ■r — < 0.14, which, in the worst scenario. 
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yields x <^ 10^ km. So, if x is not too close to nL/2, we can neglect T2 for all terrestrial 
distances. With this in mind, the formula ( p.6| ) becomes 



t{x) 



2ri 



1 + 



2ti 



(3.10) 



while linear approximation yields 



Tix) = Ti(x) . 



(3.11) 



In addition, ti can be rewritten without the term proportional to the extremely small 
coefficient uj/pq. Indeed, the numerical value of this factor in the recent experiments 
1 13, 14, ^] is less than 10~^^. Thus, with the big precision, we can write 



T"l 



ui'^ sin^ 26 sin 2ujx 
4 1 - sin^ 29 sin^ ujx 



(3.12) 



3.2 Graphic illustration 



Let us now analyze the consequences of the formula (|3.1C1|) . As we can see, the neutrino 
delay r(x) takes negative values when x E [nL, (n + ^)L], positive values when x G [{"n + 
^)L, (n + 1)L] and zeros when x = for all integers n > 0. For demonstration purposes, 
let us make a numerical example by choosing our free parameters as follows: 



Po-17 GeV, 



1.4 km. 



(3.13) 



With the known values of ArrP ^ 2.3 • 10"^ eF^, and 0.92 < sin^ 26 < 1, this yields the 
behaviour shown in Fig. ffl. The negative values of the neutrino delay time r indicate that 
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Figure 1: Oscillations of neutrino delay. The three cases are defined by sin^ 29 — 0.92 (solid line), 
sin^ 20 = 0.99 (dashed line), and sin^ 29 = 0.999 (dotted line). 



neutrino may arrive earlier than expected on the basis of the average velocity v. Having 
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in mind that u ~ 1, we expect that the neutrino velocity 



^ _ dx ^ dr(x) 
dt dx ' 



may periodically exceed the speed of light. Indeed, differentiating ( [3.10 ), we obtain 



v(x) = 1 + 4 



1 + 



1 1 - (1 + cos2 2(9) sin^ 



LOX 



1 + 



sin 26 sin 2ujx 



As expected, the neutrino velocity v{x) is a periodic function of x, with the period L. Its 
behavior is shown in Fig. ^ As we can see, most of the time it exceeds the speed of light. 
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Figure 2: Oscillations of neutrino velocity. The two cases are defined by sin^ 26 ~ 0.92 (solid line), 
and sin^ 29 = 0.95 (dashed line). 

In particular, v{x) > 1 for all x <^ L. We see that averaging v{x) over small distances 
results in superluminal effective speed. The effective speed Vefj is defined as 



^eff 



xit) 



1 



TiX 



(3.14) 



As opposed to v, that can not be directly measured, the effective speed Veff can. In fact, 
it is Veff that has been measured in the experiments ^ Using ( 3.10| ) in ( 3.14| ), 
we easily calculate the formula for the effective speed. Its graph is displayed in Fig. ^. As 
we can see, it is not exactly a periodic function, but nevertheless, it periodically exceeds 
the speed of light. 

Let us now calculate the extreme values of neutrino delays and velocities. To simplify 
exposition, we shall restrict our analysis to one period x £ [0,L]. Then, if the extreme 
points are denoted by x± {Tmin = ^i^-)-, Tmax = Tix+)), one finds x+ + x_ = L, and 
Tmax = —Tmin- The extreme points x+ and x_ are obtained by solving the equation 
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Figure 3: Oscillations of neutrino effective velocity. The three cases are defined by sin^ 2Q = 0.92 
(sohd line), sin^ 26* = 0.99 (dashed line), and sin^ 26* = 0.999 (dotted line). 



dr /dx = 0, whereupon 



L 

x± f« — 
2 



and the function Tmin = ^ix-)-, in the linear approximation in r/£, takes the value 

ujf sin2 20 

Tmin - - — . 

The corresponding effective velocity is given by Vf,ff{x_) = 1 — Tmin/x-- To illustrate this, 
let us calculate the minimal value that t{x) can possibly have for a given 9, and display the 
corresponding x_ and ^e//- Fixing our free parameters as in the example ( 3.13| ), we obtain 
the Table [l|. As we can see, the minimal time neutrino needs to arrive at x depends on 
the mixing angle 9. So does the corresponding effective velocity Veff- The maximal effect 
is obtained if 9 is close to 7r/4. Interestingly enough, the corresponding x turns out to be 



sin^ 29 


x_ 






0.92 


7, 400 km 


— 1 ns 


1 + 1.48 • 


10" 


-7 


0.99 


8, 400 km 


—3 ns 


1 + 4.22 • 


10" 


-7 


0.999 


8, 820 km 


-10 ns 


1 + 13.5 • 


10" 


-7 


0.9999 


8, 940 km 


-32 ns 


1 + 42.2 • 


10" 


-7 



Table 1: 6 dependence of the minimal neutrino delay 

close to the value x = L/2 for all the displayed 6. Unfortunately, the probability density 
Pf^ in X = L/2 goes to zero when 9 — t- 7r/4. Indeed, it is seen from ( |2.10 ) that complete 
disappearance of muon neutrino is possible only if 9 = 7r/4. Whenever this happens, a 
pure tau neutrino appears instead. 
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Similarly, the maximal value of the velocity v{x) is obtained by solving the equation 
dv/dx = 0. It is shown that Vmax increases with 6, approaching its absolute maximum 
Vmax = 1 + 1 when 9 — >• 7r/4. At the same time, x — )■ L/2. Sadly, the probability to 
detect muon neutrino in the point where it reaches its maximal velocity is zero. In other 
words, while accelerating towards Vmax — 1 + |, the muon neutrino gradually disappears. 

As for the maximum of the effective velocity fe//, it is obtained as the solution of the 
equation dv^ff/dx = 0. Again, the maximum of Vf>ff is a function of 9, which approaches 
its absolute maximum {veff)max = 1 + 2£/L in the limit 6 — )• 7r/4. At the same time, 
x — )• L/2, and the probability density goes to zero. 

Finally, let us mention that the case = 7r/4 is completely different from the case 
9 — )■ 7r/4. Indeed, when 9 = 7r/4, the velocity v{x) always exceeds the speed of light. Its 
behaviour is shown in Fig. 0. As we can see, it reaches its maximal value Vmax — 2 in 




Figure 4: Oscillations of neutrino velocity for 6 = 7r/4. The three cases are defined by uil 
(solid line), — 0.5 (dashed line), and uj£ = 0.2 (dotted line). 



0.7 



x = L/2, independently of the neutrino energy. Unfortunately, the probability density to 
detect Vmax = 2 turns out to be zero. Again, this is because the muon neutrino turns into 
tau neutrino in x = L/2. 



4. The shape of the wave packet 



Let us analyze how our results depend on the shape of the wave packet. It is useful to 
introduce some parameter 7 > 0, so that variation of 7 from zero to infinity changes the 
shape of the wave packet from almost rectangular to very sharp. Let us generalize ( |3.4D by 
choosing 



p{t) oc 



< r < 



otherwise . 



(4.1) 
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To check if the corresponding momentum distribution is sharp around p = Po, a-s required 
by our introductory assumptions, we calculate the Fourier transform 

a{p) = [ dTp{T) e*(P-Po)^ . (4.2) 



'2tt 

The resulting expression has the form 

a{po+p)^^-^, (4.3) 

where J^^i{p) are ordinary Bessel functions of the order 7 + ^. Note that a{p) is real and 

finite function of its argument. For integer values of 7, it can be expressed in terms of 
elementary functions. For example, the expression for 7 = reads 

, , sin^p ,^ ^, 

a{po +p)oc — — , (4.4) 



and the expression for 7=1 has the form 



a{po+p) (x j^il-ipcotep). (4.5) 



The necessary and sufficient condition for both amplitudes to be well localized around 
p = Po can be expressed in terms of the dimensionless quantity ipo. It reads 

£po » 1 • (4.6) 



Indeed, in the limit ipQ — )• 00, the distributions (4.4) and (4.5) take the form a{p) oc 
6{p — Po). As an illustration, the data from the recent experiments |13, 14, 15| yield 
ipo > 10^^ » 1, telling us that the condition ( [4. 6]) is indeed satisfied. In conclusion, the 
momentum distributions for the whole interval < 7 < 1 are sharply localized. 

Now that we are convinced that the approximation we have chosen to work with holds 



true, we proceed to solve the equation ( |3.lD for the set of parameter dependent amplitudes 



(4.1). Neglecting terms proportional to the small quantity T2, the equation (3^) turns into 
the quadratic equation 

TlT^ +-ffT-Tie^ = 0, 

which differs from the corresponding equation of the case 7 = 1 by the simple replacement 
Ti — >• T1/7. Having this in mind, we easily obtain 




In the lowest approximation in r/i, we come to 



r(x) = — . (4.8) 
7 
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The requirement that (4^) has a dominant role in ( [4.7D puts a restriction on 7. For the 
experiments |13, 14, a sufficient condition to ensure the validity of ( [4.8D is 



7 > 10 



-8 



Finally, let us say something about how neutrino velocity depends on 7. We have 
seen in the preceding section that Vmax is a discontinuous function in 6 = n / A. The same 
happens in the case 7 7^ 1. Indeed, the neutrino velocity has two absolute maximums, 



27 

1 H when 

3 

= 1+7 when 



TT 

' 4 

vr 



which both take place in x = L/2. In this point, however, the probability density goes 
to zero. Thus, we can say that muon neutrino, which approaches the point of absolutely 
maximal velocity, gradually disappears. This is because it turns into a tau neutrino. 




Figure 5: Shapes of the neutrino wave packet. The three cases are defined by 7 = 1 (solid line), 
7 = 0.1 (dashed line), and 7 — 0.01 (dotted line). 



In conclusion, as 7 approaches zero, the wave packet is more sharply localized in 
momentum space, and has more rectangular shape in coordinate space. (See the example 
in Fig. |5|.) This makes the neutrino delay time |r| increase. As a consequence, the 
superluminal effect becomes easier to detect. Still, as the known experiments lack the 
parameter that determines the shape of the neutrino wave packet, the comparison with 
measurements remains an uneasy problem. 

At the end, let us note that the above results refer to the detection of muon neutrino. 
The generalization to the case of tau neutrinos is straightforward. In the next section, we 
shall analyze the evolution of tau component of initially pure muon neutrino. 
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5. Velocity oscillations of tau neutrinos 



In the preceding section, we have considered the dynamics of initially pure muon neutrino, 
and calculated the probability to detect muon neutrino at later times. In a similar way, we 
can derive the probability density to detect tau neutrino. It is defined by 

Pr{x,t) = \{iyr{x)\iy^{t))f , 

where |z^r(2;)) = sin 6* + cos 6 \ 2)\x) is the eigenstate of the position operator for tau 
neutrino. The direct calculation yields 

Pr{x,t) = p^{vT)s\n^2e s\n^ujt. (5.1) 

As we can see, no correction linear in Aw appears in ( |5.lD . It is checked that 

j {P^ + Pr) dx = j p'^{x)dx = 1 , 

as it should be in the approximation we work with. 

If the neutrino detector detects tau neutrinos instead of muon neutrinos, the corre- 
sponding equation dP-r/dt = takes the form 

^ + ^-3(x) = 0, (5.2) 

where 

T2,{x) = cotujx . (5-3) 



As we can see, the dependence on the mixing angle 6 does not appear in ( |5.3| ). Choosing 
the same form for p as in ( [4.1[) , we obtain the expression of the same form as ( |4.7] ), with 
the only difference that t^{x) is used instead of ti{x). The approximation linear in t/£ 
then yields 

r.(x) = ^. (5.4) 
7 

If we restrict our analysis to distances which are not too close to (n + |)L, and the 
mixing angle 6 is close to 7r/4, the corresponding formula for ti takes the simplified form 
n w — ^ tanwx. It then yields 



Ti[x) T-i[x) = const. , (5.5) 



or equivalently. 



a;2 



T,{x) Trix) « [f - r2(x)] [f - r'Ax)] . (5.6) 

As we can see, the arrival times and Tr for muon and tau neutrinos have opposite signs. 
Thus, when one of them has superluminal speed, the speed of the other is subluminal. The 
relation (|5.6| ) nicely illustrates connection between oscillations of muon and tau neutrinos. 
When one fiavor has minimal deviation (r^(x) — >• 0) the other has maximal one (t,-(x) — >• i). 
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(Note, however, that our results are derived under the assumption that r is not too close 
either to I or to zero.) 

Finally, let us point out another interesting relation between t\(x) and r-iix): 

Ti(x + L/2) «2t3(x). 

With this, the calculations related to tau neutrino are greatly simplified. In particular, the 
linear approximation in r/^ yields the relation t^{x + L/2) = 2Tr{x), showing that the 
double delay time of tau neutrino is obtained by the simple shifting x — )• x + L/2 in the 
expression for the delay time of muon neutrino. 



6. Summary and discussion 

In this paper, we have demonstrated that the flight of a free neutrino is characterized not 
only by the well known flavor oscillations, but also by the oscillations of neutrino velocity. 
This has been done by considering the free evolution of initially pure muon neutrino. First, 
we calculated the probability density for detecting muon neutrino by the detector placed in 
a fixed position x. Such a probability distribution was a function of time, and its maximum 
was identified with the moment neutrino arrived at x. This way, we obtained the formula 
for the evaluation of time the neutrino needed to fly across the distance x. 

It should be noted that our formula differs from similar formulae found in literature 
(see, for example, 11, [l^). The reason for this is the difference in our definitions of the 



neutrino position. While we identify it with the maximum of its probability distribution. 



the authors of 11, 12j use the average of neutrino position operator. As a consequence, 
our formula carries additional dependence on the size and shape of the neutrino wave 
packet. The neutrino delay, as compared to the arrival time of the photon, was found to be 
a periodic function of x, with the period which did not depend on the size and shape of the 
neutrino wave packet. The neutrino velocity, on the other hand, turned out to drastically 
depend on the size and shape of the wave packet: the more rectangular the shape of the 
packet was, the higher velocity neutrino could reach. As it turned out, the neutrino velocity 
periodically exceeded the speed of light. We were able to derive the maximal velocity the 
neutrino could possibly achieve during its flight. As it turned out, the probability to detect 
maximal velocity was rather small, approaching zero when 9 — )• tt/4. 

The fact that our formula depends on the size and shape of the neutrino wave packet 
makes the comparison with the experiment very difficult. This is because the known 
experiments do not provide the information on these two parameters. Nevertheless, we 
shall try to test our formula by comparing its predictions with three recent experiments 
13, 14, 15|. To this end, we shall estimate the size of the neutrino wave packet using the 



results found in literature. In particular, a useful study of the properties of accelerator 



neutrinos can be found in Ref. |T^. Using the size of the wave packet as found in [16|, 
and almost Gaussian shape as defined by 7 = 1, we make the Table |2[ As we can 
see, our theoretical predictions do not contradict any of the three experiments. Even if 
the neutrino wave packet is taken 10^ times longer, the agreement with the experiment 
is not compromised. The same holds when it comes to the wave packet shape. Indeed, 
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Experiment 


Po 


2£ 


X 


Vexp 1 


Veff - 1 


MINOS 
ICARUS 
OPERA 


3 GeV 
17 GeV 
17 GeV 


7.7 cm 
0.67 cm 
0.67 cm 


734 km 
730 km 
730 km 


(5.1 ±7.5) • 10-5 
(0.1 ±5.7) • 10-6 
(2.7 ±6.5) • 10-6 


0.9- 10-15 
1.5 • 10-19 
1.5 • 10-19 



Table 2: Comparison witii experiments 



the change of 7 in the ahowed interval 7 ^ 10-^ leaves the theoretical values within the 
experimental error. In fact, it is the big experimental error that basically ensures this 
agreement. For a real test of our theoretical predictions, more efficient measurements are 
needed. In particular, our equations suggest how neutrino free parameters should be chosen 
to maximize the superluminal effect. 

Let us say something about energy dependence of the neutrino velocity. As seen from 
our formulae, the energy dependence of v^ff — 1 has oscillating character. To simplify the 
analysis, we shall restrict to the region x <^ L, which is achieved by using short range high- 
energy neutrinos. Upon this, the measured quantity Veff — 1 becomes proportional to £'^/pq, 
which reduces to i'^/E^ in the ultrarelativistic limit. At the same time, the wave packet 
size i is also energy dependent. Indeed, it has been shown in [16| that, depending on the 
experimental details, £ is proportional to either 1/Eq or 1/Eq. Thus, in the ultrarelativistic 
limit, 

1 1 

^'^^ ^ eI ^^^-^ ^ eI' 

This is a strong energy dependence, but still undetectable by the ICARUS and OPERA 
experiments. (MINOS experiment is excluded from the analysis because it violates the 
requirement x <^ L.) Indeed, it is seen from Table |^ that the values of Ue// — 1 are IQi^ 
times smaller than Vexp — 1, which implies that 100 times lower energy still yields the result 
that agrees with the experiments. Higher energies, on the other hand, diminish the value 
of Vf.ff — 1, and therefore, fit the experiments even better. 

Finally, let us comment on the result of Ref. |17] stating that superluminal neutrinos 
should rapidly loose energy during their flight. This result has been derived with the 
assumption that neutrino speed does not change along neutrino trajectory. However, we 
have shown that this is not the case. In general, the neutrino velocity has oscillating 
character, but we shall simplify the analysis by adopting the condition x <^ L. In this 
regime, the results of the preceding sections yield 

d f X \ 1 f /\m^t ^ 

2a— — 7T , a = — — — 



dxKE'^J' 27 V 

where 6 = v'^ — 1. Using this expression in the formula for the rate at which superluminal 
neutrinos loose their energy |17], 

dE .,0 25 Gl 



dx ' 448 1927r3 ' 

we obtain a complicated, higher order differential equation. To simplify the calculations, 
we shall assume that superluminal neutrinos loose their energy slowly (\dE/dx\ <C E/x). 
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Then, we obtain a simple approximative solution 



E^Eq-Et, Et = Sko^ X , 

telling us that the rate superluminal neutrinos loose their energy is linear in x. Notice, 
however, that the value of the constant na^ is extremely small in the present experiments. 
For the numerical illustration, we shall take data from Table ^. With the known value of 
the Fermi coupling constant Gp ~ 1.17- 10^^ GeV~^, and using the almost Gaussian wave 
packet as defined by 7 = 1, we find 

Et < 1.1 • 10~^^ GeV. 



Having in mind that the initial neutrino energy in the considered experiments [fL3| , 14, |15[ 
is Eq > 3 GeV, we see that the loss of energy during the flight of superluminal neutrinos is 
negligible. Our initial assumption \dE/dx\ <C E/x is thereby justified, and the validity of 
our conclusion is confirmed. We can still change the wave packet shape, but this can not 
significantly modify our conclusions. Indeed, the lowest allowed value of 7 is of the order 
10"'^, which leads to the value Et ~ 10~^^ GeV. This is still negligible with respect to 
the initial neutrino energy (Et <^ Eq) in all terestrial experiments. Summarized, we have 
proven the existence of energy conserving superluminal free neutrinos. 

In conclusion, if neutrinos have different masses, the oscillations of neutrino velocity 
necessarily exist. In particular, the neutrino velocity periodically exceeds the speed of 
light. The significance of this result is threefold. First, it shows that superluminal speed 
can be achieved without violation of special relativity. Second, our equations suggest 
how neutrino parameters should be chosen to maximize the superluminal effect in new 
experiments. Finally, our formula offers an independent way to determine neutrino mixing 
angles. 
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